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ABSTRACT 

We have analyzed the clustering of C IV absorption line systems in an exten- 
sive new catalog of heavy element QSO absorbers. The catalog permits explo- 
ration of clustering over a large range in both scale (from about 1 to over 300 
hr x Mpc) and redshift (z from 1.2 to 4.5). We find significant evidence (5.0 a; 
Q = 2.9 x 10~ 7 ) that C IV absorbers are clustered on comoving scales of 100 
hr x Mpc (g = 0.5) and less — similar to the size of voids and walls found in 
galaxy redshift surveys of the local universe (z < 0.2) — with a mean correlation 
function (£ aa ) = 0.42 ± 0.10 over these scales. We find, on these scales, that the 
mean correlation function at low ((z)\ ow = 1.7), medium ((z) me d = 2.4), and high 
redshift {(z) high = 3.0) is (£ aa ) = 0.40 ±0.17, 0.32 ±0.14, and 0.72 ±0.25, respec- 
tively. Thus, the superclustering is present even at high redshift; furthermore, 
it does not appear that the superclustering scale, in comoving coordinates, has 
changed significantly since then. We find 7 QSOs with rich groups of absorbers 
(potential superclusters) that account for a significant portion of the clustering 
signal, with 2 at redshift z ~ 2.8. We find that the superclustering is just as 
evident if we take g = 0.1 instead of 0.5; however, the inferred scale of clustering 
is then 240 hr x Mpc, which is larger than the largest scales of clustering known 
at present. This discrepancy may be indicative of a larger value of g , and hence 
f2 . The evolution of the correlation function on 50 Yc 1 Mpc scales is consistent 
with that expected in cosmologies with density parameter ranging from Qq — 
0.1 to 1. Finally, we find no evidence for clustering on scales greater than 100 
h- 1 Mpc (g = 0.5) or 240 hr 1 Mpc (g = 0.1). 

Subject headings: catalogs — cosmology: observations — large-scale structure 
of universe — quasars: absorption lines 
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1. Introduction 

It has been recognized for some time now that QSO absorption line systems are partic- 
ularly effective probes of large-scale structure in the universe (see, e.g., |Shaver fc Robertson 



|1983| , |Crotts, Melott, fc York 1985|) . This is because the absorbers trace matter lying on the 



QSO line of sight, which can extend over a sizable redshift interval out to high redshifts. 
Thus, the absorbers trace both the large-scale structure (on scales out to hundreds of Mpc) 
and its evolution in time, since the clustering pattern can be examined as a function of 
redshift out to z ~ 4. The evolution of large-scale structure is of great interest, since, in the 
gravitational instability picture, it depends sensitively on the mean mass density Qq ( Peebles 
19801 , 1993). 



We are not concerned here with the relationship between the absorbers and galaxy 
haloes, or with small-scale structure. However, in another paper ( IQuashnock fc Vanden Berk] 



1996| ) we discuss the small-scale clustering of absorbers and relate it to galaxy clustering. 



In this Letter, we present results of an analysis of line-of-sight correlations of C IV 
absorption line systems, using a new and extensive catalog of absorbers ( |Vanden Berk et 



|al. 1996a]) . This catalog contains data on all QSO heavy-element absorption lines in the 
literature, complete up to December 1995, with some additional entries since then. It is 
an updated version of the catalog of York et al. (1991), but is more than twice the size, 
with over 2200 absorbers listed over 500 QSOs, and is the largest sample of heavy-element 
absorbers compiled to date. More details can be found in the earlier version of the catalog, 
as well as in our recent paper flVanden Berk et al. 1996T5| ) in which we find that the number 



of absorbers is correlated with the intrinsic brightness of the QSO, suggesting that QSOs 
are lensed by matter associated with the absorbers. A preliminary study suggests that the 
absorber correlation function is not strongly dependent on the intrinsic brightness of the 
QSO. We describe the correlation analysis in §2, present our results in §3, and discuss the 
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implications of these in §4. 



2. Correlation Analysis of C IV Absorbers 

Here we briefly describe our procedure for calculating the line-of-sight correlation func- 
tion, £ aa . Unless otherwise noted, we take go = 0.5 and A = 0. We follow the usual convention 
and take h to be the Hubble constant in units of 100 km s _1 Mpc -1 . 

To produce a more uniform set of absorbers from the inhomogeneous catalog, we have 
applied the following selection criteria to the QSO spectra and absorbers. The Lya forest 
region of each spectrum was excluded because identification of heavy-element lines there 
is problematic. This means we examine absorbers within about 60,000 km s _1 , or about 
400 h^ 1 Mpc, of the QSO. The so-called associated region within 5000 km s _1 of each QSO 
emission redshift has also been excluded, since the number density of absorbers there may not 
be representative of that for absorbers farther removed in redshift from their QSOs flAldcroft 



Bechtold, fc Elvis 1994j) . Thus a typical redshift range for a line of sight is Az = 0.4. All 



of the selected absorption systems must at least have an identified C IV doublet, and if no 
other line was identified, the equivalent width ratio of the 1548 A component to the 1551 A 
component must be > 1 within the listed measurement errors. In addition, the equivalent 
width of the 1548 A component must have been detected at more than the 5 a level, and have 
a rest value of at least 0.1 A. In order to avoid aliasing of power ([Tytler, Sandoval, fc Fan 



1993| ) on the large scales of interest because of the very pronounced peak in the correlation 



function on small scales QSargent, Boksenberg, fc Steidel 19881 , IQuashnock fc Vanden Berk 



1996 ), we have combined all absorbers lying within 3.5 comoving h~ l Mpc of each other 
(~ 600 km s _1 ) into a single system with redshift equal to the average redshift of the 
components and with equivalent width equal to that of the strongest component. Because of 
this, we have included only spectra with resolutions that allow absorber comoving separations 
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< 3.5 h Mpc to be resolved. 

These selection criteria leave a total of 360 C IV absorption systems, with redshifts from 
1.2 to 4.5, scattered among 373 QSO lines of sight. The distribution in the total number 
of absorbers over the line of sight is consistent with that from a uniform distribution. For 
number of absorbers per line of sight through 7, the observed versus expected number of 
lines of sight are (207, 190.34), (67, 82.52), (49, 52.23), (22, 28.24), (16, 12.92), (9, 4.61), (1, 
1.60), and (2, 0.39). A KS-test shows that the two distributions are consistent with each 
other (Q = 0.45). This shows that the average distribution of absorbers on 400 h^ 1 Mpc 
scales is uniform. 

To study the distribution on smaller scales, we use the C IV line-of-sight correlation 
function, £ aa , which is calculated by binning the comoving separation, Ar, of absorber pairs 
and comparing the number of real pairs, iV n in a bin to the average number of pairs, N s , 
generated from 100 Monte Carlo simulations; namely, £ aa (Ar) = (N r /N s ) — 1. The simulated 
absorber samples are generated by randomly redistributing all of the real absorbers among 
the QSO spectra, subject to the condition that the absorbers could have been detected in 
both redshift and equivalent width in their randomly selected QSO spectra. This technique 
is possible because the catalog of absorbers lists the wavelength limits and equivalent width 
limits for each of the listed QSO spectra. This technique has the virtue that the redshift 
and equivalent width number densities of all the simulated data sets are exactly the same 
as for the real data set, so that the real and simulated data sets have the same evolution 
of physical properties as a function of redshift. The 1 a region of scatter around the null 
hypothesis of no clustering (£ aa = 0) is given by the standard deviation of £ aa in the 100 
random simulations, whereas the 1 a error in the estimator of £ aa in each Ar bin is N s 
( [Peebles 1980| ). When clustering is present, these errors will be different. 
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3. Results 

The size and extent of the absorber catalog of Vanden Berk et al. (1996a) permits 
exploration of clustering over an unprecedented range in scale (from about 1 to over 300 
h^ 1 Mpc) and redshift (z from 1.2 to 4.5). Figure 1 shows the line-of-sight correlation 
function of C IV absorbers, £ aa (Ar), as a function of absorber comoving separation, Ar, for 
the entire sample of absorbers. The results are shown for both a go — 0.5 (top panel, 25 
h^ 1 Mpc bins) and a go — 0.1 (bottom panel, 60 h -1 Mpc bins) cosmology.^] The vertical error 
bars through the data points are 1 o errors in the estimator for £ aa , which differ from the 
1 a region of scatter (dashed line, calculated by Monte Carlo using the bootstrap technique 
described above) around the no-clustering null hypothesis. 

Remarkably, there appears to be significant clustering in the first four bins of Figure 1: 
The mean correlation for those bins is (£ aa ) = 0.42 ± 0.10 (go = 0.5) or 0.39 ± 0.10 (go = 
0.1). Here the errors quoted are the la errors in the estimate of the correlation function 
(corresponding to the error bars in Figure 1). In assessing the significance of these results, 
however, the appropriate measure of departure from uniformity is the 1 a region of scatter 
around the null hypothesis (corresponding to the dashed lines in Figure 1), which equals 
0.085 (go = 0.5) or 0.082 (go = 0.1). Thus, the positive correlation seen in the first four bins 
of Figure 1 has a significance of 5.0 a (Q = 2.9 x 10 -7 ), if g = 0.5, and 4.8 a (Q = 9.2 x 10 -7 ) 
if go = 0.1. Therefore, there is significant evidence of clustering of matter traced by C IV 
absorbers on scales up to 100 hr x Mpc (g = 0.5) or 240 hr x Mpc (g = 0.1). 

There is no evidence from Figure 1 for clustering on comoving scales greater than these. 
The mean correlation of the last six bins in Figure 1 is (£ aa ) = 0.01 ± 0.08 (go = 0.5 or 0.1), 

2 Larger bins are required for go = 0.1 because, at high redshift, a larger comoving sepa- 
ration Ar arises from a fixed redshift interval Az. 
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which is consistent with zero. 

We have investigated the evolution of the superclustering by dividing the absorber sam- 
ple into three approximately equal redshift sub-samples. Figure 2 shows £ aa (Ar), as a 
function of absorber comoving separation, Ar (q = 0.5), for low (1.2 < z < 2.0, top panel), 
medium (2.0 < z < 2.8, middle panel), and high (2.8 < z < 4.5, bottom panel) redshift ab- 
sorber sub-samples. Larger bins for the comoving separation were used in Figure 2 because 
of the smaller number of absorbers in each sub-sample. The error bars are as in Figure 1. 

The correlation function estimator is positive in the first two bins (corresponding to 
scales of 100 hr l Mpc and less, like the first four bins in Figure 1, top panel) in all three 
panels of Figure 2: The mean correlation (q = 0.5) is (£ aa ) = 0.40 ± 0.17 (2.7 a) for the 
low redshift sub-sample, (£ aa ) = 0.32 ±0.14 (2.4 a) for the medium redshift sub-sample, and 
(£aa) = 0.72 ± 0.25 (3.4 a) for the high redshift sub-sample. The values for q = 0.1 are 
essentially the same as these. 

Thus, the significant superclustering seen in Figure 1 is present in all three redshift 
sub-samples in Figure 2, so that the superclustering is present even at redshift z ^ 3. 
Furthermore, it does not appear that the superclustering scale, in comoving coordinates, has 
changed significantly since then. 

4. Discussion 

We have found evidence for the existence of large-scale superclustering that has existed 
since redshift -2^3 and has a comoving scale ~ 100 hr l Mpc (q = 0.5) or 240 h~ l Mpc 
(go = 0.1) that has not changed much since then. We have examined the clustering signal 
more closely and find that a large portion comes from 7 QSO lines of sight that have groups 
of 4 or more C IV absorbers within a 100 hr x Mpc interval (q = 0.5). (From Monte Carlo 
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simulations, we expect only 2.7 ±1.5 QSOs with such groups.) When these are removed, the 
mean correlation function on scales of 100 h^ 1 Mpc and less (go — 0.5) is (£ aa ) = 0.16 ±0.14. 
Table 1 lists these 7 QSOs, along with the properties of their groups of C IV absorbers. 

Several of the groups in Table 1 have been previously identified as potential superclusters 
( |Romani, Filippenko, fc Steidel 1991| , [Aragon-Salamanca et al. 1994Q ; indeed, superclustering 



has been suggested because of an unusual concentration of absorbers, within several tens of 
megaparsecs, in the lines of sight of 1037-2704 and 1038-2712 flJakobsen et al. 1986| , |S argent 



fe Steidel 1987| , [Robertson 1987| , Pinshaw fc Impey 1996| , |Lespine fc Petitjean 1996j ), and 



0237-233 flBoisse 1987j , [Heisler, Hogan, fc WhiteT^ , [Koltz et al. Tggg). However, there 



was concern that any inferred superclustering might be overestimated because of aliasing 
of small-scale power. Like Dinshaw & Impey (1996), we find that superclustering is not 
caused by aliasing, since it remains after combining all absorbers lying within 3.5 h" 1 Mpc 
of each other into a single system. (For example, we count the 11 absorbers of 0237-233 
in a 100 h~ l Mpc interval as being 4 independent systems.) We have found two potential 
superclusters in the spectra of QSOs 2126-158 and 2359±068, at redshift z ~ 2.8. 

The superclustering is indicative of generic large-scale clustering in the universe, out to 
high redshift z ^ 3, on a scale frozen in comoving coordinates that is — if go — 0.5 — similar 
to the size of the voids and walls in galaxy redshift surveys of the local universe (z < 0.2), 
such as the CfA survey ( IGeller fc Huchra 1989Q and others ( [Kirshner et al. 19811 , |Ch mcarim 



Oiovanelli, fc Haynes 1983] , |da Uosta et al. 1994| , |Landy et al. 1996] ) . It also appears consistent 



with the general finding of Broadhurst et al. (1990) that galaxies are clustered on very large 
scales, although we have not confirmed that there is quasi-periodic clustering with power 
peaked at 128 h' 1 Mpc. It may also have been found in quasars ( [Kravtsov 1996|) . 

Our estimate of the superclustering scale increases to 240 hr 1 Mpc if go = 0.1 (see 
Figure 1), which is larger than the largest scales of clustering known at present. If the 
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structures traced by C IV absorbers are of the same nature as those seen locally in galaxy 
redshift surveys, the superclustering scale should have a value closer to 100 h^ 1 Mpc . This 
may be indicative of a larger value of qo, and hence £Iq. 

In Figure 3, we show the mean correlation function, (£ aa ), on scales of 50 h~ x Mpc and 
less, as a function of redshift z, for the entire data set of absorbers (solid symbols), and for 
strong Mg II absorbers ( (Steidel fc Sargent 1992| , open symbol), with go — 0.5. Also shown is 



its expected dependence on redshift, calculated from the growth factor ( Peebles 1980| , 1993), 
normalized to the sample mean correlation of 0.48 at the sample median redshift of z = 2.31, 
for Qq — 1 (solid line) and Qq = 0.1 (dashed line). While the evolution of the correlation 
function on 50 h~ l Mpc scales is consistent with that expected in cosmologies with density 
parameter ranging from Q = 0.1 to 1, the Sloan Digital Sky Survey should give better 
estimates for (£ aa ) (see the typical error bar at the upper right of Figure 3) and will be able 
to constrain Qq to 20%. 

It is possible that structures of size ~ 100 h^ 1 Mpc may have arisen through the 
Zel'dovich "pancake" scenario ( [Zerdovich 1970| , Melott fc Shandarin 1990| , [Peebles 1993| ). 



and that this particular scale has deep significance because it is of order the horizon size at 
decoupling, and is the scale on which an "acoustic peak" is expected in the initial power 
spectrum of many cosmologies ( Peebles 1993 . [Landy et al. 19961) . The lack of clustering on 



scales greater than 100 h 1 Mpc has implications on the peak and turn-over scale of the 
power spectrum ( |Quashnock fc Vanden Berk 1996| ). 



We acknowledge helpful discussions with Don Lamb, Paul Ricker, Bob Rosner, and Alex 
Szalay, and statistical comments by Carlo Graziani. Damian Bruni and Chris Mallouris have 
greatly assisted in compiling the new catalog. JMQ is supported by the Compton Fellowship 
- NASA grant GDP93-08. DEVB was supported in part by the Adler Fellowship at the 
University of Chicago, and by NASA Space Telescope grant GO-06007.01-94A. 
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Table 1 


OSOs with 4 or 


more C IV absorbers in a 


100 h~ l Mnc 


infprvpu ( Hc\ — D Fi\ 


QSO 


N 

1 v exp 


iVabs 


(^abs) 


A^abs 


Ar (h- 1 Mpc) 


0149+336 


0.57 


4 


2.1583 


0.1593 


85.8 


0237-233 


0.30 


4 


1.6348 


0.0775 


54.4 


0958+551 


0.25 


4 


1.3207 


0.1045 


88.5 


1037-270 


0.90 


5 


2.0133 


0.1703 


98.5 


1038-272 


0.79 


4 


1.9283 


0.1642 


98.2 


2126-158 


0.47 


5 


2.7265 


0.1814 


75.6 


2126-158* 


0.74 


5 


2.7804 


0.2280 


92.6 


2359+068 


0.51 


6 


2.8303 


0.2125 


84.8 


2359+068* 


0.64 


6 


2.8730 


0.2394 


94.5 



*Part of this group forms another group (also listed) with an additional absorber within 
the same QSO, but the extent of the two groups taken together exceeds 100 h~ l Mpc. iV ex p 
is the number of expected absorbers in the interval Az a b s - 
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Fig. 1. — Line-of-sight correlation function of C IV absorbers, £ aa (Ar), as a function of 
absorber comoving separation, Ar, for the entire sample of absorbers, in a go = 0.5 (top 
panel, 25 h~ l Mpc bins) and a go = 0.1 (bottom panel, 60 h~ l Mpc bins) cosmology. The 
vertical error bars through the data points are 1 a errors in the estimator for £ aa , which 
differ from the 1 a region of scatter (dashed line) around the null hypothesis of no clustering 

(U = 0). 

Fig. 2. — Line-of-sight correlation function of C IV absorbers, £ aa (Ar), as a function of 
absorber comoving separation, Ar (go = 0.5), for low (1.2 < z < 2.0, top panel), medium 
(2.0 < z < 2.8, middle panel), and high (2.8 < z < 4.5, bottom panel) redshift absorber 
sub-samples . The vertical error bars through the data points are 1 a errors in the estimator 
for £ aa . 

Fig. 3. — Mean correlation function, (£ aa ), on scales of 50 h~ x Mpc and less, as a function 
of redshift z, for the entire data set of C IV absorbers (solid symbols), and for strong Mg II 
absorbers flSteidel fc Sargent 1992| , open symbol), with g = 0.5. Also shown is its expected 



dependence on redshift, normalized to the sample mean of 0.48 at the median sample redshift 
of z — 2.31, for = 1 (solid line) and Qq = 0.1 (dashed line). The expected error bar from 
the Sloan Digital Sky Survey is shown in the upper right. 
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